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Abstract  -  In  the  defence  domain ,  weapons  allocation  is  de¬ 
fined  to  be  the  reactive  assignment  of  weapon  systems  to  engage 
or  counter  identified  threats.  From  a  command  perspective,  this 
refers  to  the  allocation  of  friendly  and  coalition  force  elements 
(eg  fighter  aircraft,  frigates  etc )  to  engage  or  interdict  adver¬ 
saries  which  are  posing  threats,  not  only  to  themselves,  but  also 
to  defended  areas  and  high-value  assets.  In  an  earlier  work, 
a  conceptual  rule-based  approach  to  weapons  allocation  in  the 
air  domain  was  outlined  in  terms  of  so-called  critical  and  sort¬ 
ing  parameters,  which  may  be  used  to  determine  the  capability 
of  each  friendly  airborne  interceptor  to  engage  or  counter  each 
threat,  and  to  rank  the  candidate  interceptor-threat  pairings  re¬ 
spectively.  An  issue  of  relevance  to  the  evaluation  of  the  (in¬ 
terdependent)  parameters  of  fuel  sufficiency,  egress  safety  and 
time-to-intercept  is  how  to  determine  the  shortest  path  from  a 
given  interceptor  to  a  static  or  dynamic  threat  which  avoids  pro¬ 
hibited  areas  such  as  missile  engagement  zones,  neutral  and  en¬ 
emy  territories  and  other  exclusion  zones.  However,  in  general 
finding  the  shortest  path  is  a  non-trivial  exercise  and  so  deter¬ 
mining  suboptimal  paths  through  the  prohibited  areas  is  often 
necessary.  In  the  current  paper,  the  problem  is  investigated  from 
both  perspectives.  In  particular,  a  technique  developed  for,  and 
applied  to,  the  field  of  robotics  for  finding  the  shortest  path  from 
a  source  to  a  fixed  destination  through  a  flat  earth  environment 
littered  with  obstacles  is  adapted  to  solve  the  shortest  path  prob¬ 
lem  for  a  spherical  earth  geometry.  This  is  then  used  as  the  basis 
for  determining  efficient  paths  from  an  interceptor  to  engage  or 
counter  a  moving  threat. 

Keywords:  Weapons  allocation,  threat  evaluation,  air  defence, 
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1  Introduction 

In  the  defence  domain,  weapons  allocation  refers  to  the  re¬ 
active  assignment  of  weapon  systems  to  engage  or  counter 
identified  threats.  From  a  single  platform  perspective,  this 
typically  refers  to  the  allocation  of  the  platform’s  arma¬ 
ments  (eg  missiles,  bombs,  torpedoes  etc)  or  countermea¬ 
sures  (eg  chaff,  flares,  decoys  etc)  to  engage  or  evade  the 
various  threats  facing  it  (see  for  example  [1]).  From  a  com¬ 
mand  perspective,  however,  this  refers  more  to  the  alloca¬ 
tion  of  friendly  and  coalition  force  elements  (eg  fighter  air- 
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craft,  frigates  etc)  to  engage  or  interdict  adversaries  which 
are  posing  threats,  not  only  to  themselves  and  other  force 
elements,  but  also  to  defended  areas  and  high-value  ground 
and  surface-based  assets.  It  is  the  latter  perspective  which 
is  relevant  to  the  current  paper. 

In  an  earlier  work  [2],  a  conceptual  rule-based  approach 
to  weapons  allocation  in  the  air  domain  has  been  outlined 
in  terms  of  so-called  critical  and  sorting  parameters,  which 
may  be  used  to  determine  the  capability  of  each  friendly 
airborne  interceptor  to  engage  or  counter  each  threat,  and  to 
rank  the  candidate  interceptor-threat  pairings  respectively. 
It  lists  the  critical  parameters  as: 

•  Threat  priority  -  This  determines  the  order  in  which 
threats  are  dealt  with  from  the  most  to  the  least  threat¬ 
ening  (if  they  are  to  be  handled  sequentially); 

•  Interceptor  availability  -  This  indicates  whether  the  in¬ 
terceptor  is  available  for  allocation  to  a  threat  or  if  in¬ 
stead  it  has  already  been  assigned  a  task; 

•  Interceptor  suitability  -  This  indicates  whether  the 
platform  is  of  the  correct  category  eg  fighter  or 
bomber,  for  allocation  to  a  threat; 

•  Weapons  effectiveness  -  This  determines  whether  the 
interceptor’s  weapons  are  capable  of  defeating  the 
threat; 

•  Fuel  sufficiency  -  This  determines  if  the  intercep¬ 
tor  has  sufficient  fuel  to  reach  the  threat,  engage  or 
counter  it  and  then  return  to  base  or  force;  and 

•  Egress  safety  -  This  determines  if  the  interceptor  will 
have  an  unobstructed  path  on  its  return  to  base  or  force, 
after  having  engaged  or  countered  the  threat, 

while  the  sorting  parameters  are: 

•  Time-to-intercept1  -  This  is  an  estimate  of  the  time  re¬ 
quired  for  the  interceptor  to  come  into  range  of  the 
threat  to  engage  or  counter  it; 

•  Intercept  path  -  This  determines  whether  the  intercep¬ 
tor  may  intercept  the  threat  without  having  to  pass 
through  hostile  territory  or  other  prohibited  areas;  and 

1  Arguably,  time-to-intercept  should  also  be  used  as  a  critical  parameter, 
since  an  interceptor  should  not  be  considered  for  pairing  with  a  threat  if  it 
can’t  reach  it  before  it  carries  out  its  mission. 
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•  Intercept  angle  -  Based  on  the  required  nature  of  the 
intercept,  this  determines  if  the  intercept  geometry  is 
suitable  (eg  collision-collision  intercepts  for  engage¬ 
ments  and  pure  pursuit  intercepts  for  identification). 

An  issue  of  relevance  to  the  evaluation  of  the  (interdepen¬ 
dent)  parameters  of  fuel  sufficiency,  egress  safety  and  time- 
to-intercept  is  how  to  determine  the  shortest  path  from  a 
given  interceptor  to  a  static  or  dynamic  threat  which  avoids 
prohibited  areas  such  as  missile  engagement  zones,  neu¬ 
tral  territories  and  other  exclusion  zones.  However,  in  gen¬ 
eral,  finding  the  shortest  path  is  a  non-trivial  exercise  and  so 
determining  suboptimal  paths  through  the  prohibited  areas 
which  are  efficient  in  terms  of  the  distance  travelled  and  the 
time  to  earliest  weapon  release  is  often  necessary.  In  the 
current  paper,  the  problem  is  investigated  from  both  per¬ 
spectives. 

In  particular,  a  motion  planning  technique  developed  for, 
and  applied  to,  the  field  of  robotics  for  finding  the  shortest 
path  from  a  source  to  a  fixed  destination  through  a  flat  earth 
environment  littered  with  obstacles  is  adapted  to  solve  the 
shortest  path  problem  for  a  spherical  earth  geometry.  This 
is  then  used  as  the  basis  for  determining  the  shortest  path 
from  the  interceptor  to  engage  or  counter  a  static  threat  and 
an  efficient  path  from  an  interceptor  to  engage  or  counter  a 
moving  threat. 

The  remainder  of  the  paper  is  structured  as  follows.  In 
Section  2,  the  threat  interception  problem  is  formulated.  In 
addition,  relevant  mathematical  notation,  terminology  and 
results  are  introduced.  In  Section  3,  the  threat  interception 
problem  in  the  absence  of  prohibited  areas  is  discussed  for  a 
spherical  earth  geometry.  In  Section  4,  the  so-called  short¬ 
est  path  problem  is  discussed  in  detail.  In  Section  5,  the 
solution  to  the  shortest  path  problem  is  used  to  modify  the 
approach  to  the  threat  interception  problem  to  allow  for  the 
presence  of  prohibited  areas.  In  Section  6,  a  brief  example 
is  presented  to  illustrate  the  proposed  techniques.  Finally, 
some  concluding  remarks  are  made  on  possible  avenues  for 
further  research. 

2  Problem  Formulation  and  Mathe¬ 
matical  Details 

In  the  context  of  tactical  air  defence,  consider  a  scenario  in 
which  a  threat  evaluation  has  determined  that  a  given  threat 
T  is  to  be  intercepted.  Furthermore,  suppose  that  a  friendly 
interceptor  I  is  a  candidate  for  the  interception.  In  evalu¬ 
ating  the  pairing  of  /  to  T,  it  needs  to  be  determined  if  I 
is  capable  of  intercepting  T,  that  is  bringing  T  within  its 
weapons  envelope,  and  if  so,  what  the  minimum  time  to  the 
intercept  (also  know  as  the  time  to  earliest  weapon  release) 
is  and  what  path  I  should  follow  to  achieve  the  intercep¬ 
tion  within  this  minimum  time.  Allowance  also  needs  to 
be  made  for  the  fact  that  I  is  to  avoid  flying  through  the 
airspace  of  nominated  prohibited  areas  such  as  missile  en¬ 
gagement  zones,  neutral  and  enemy  territories  and  other  ex¬ 
clusion  zones  for  example.  A  depiction  of  such  a  scenario 
with  respect  to  a  spherical  earth  geometry  appears  in  Fig.  1, 
featuring  the  interceptor  on  the  left  in  blue,  the  threat  (and 
its  initial  heading)  on  the  right  in  red,  and  the  prohibited 
areas  between  them  in  magenta. 


Figure  1:  A  Threat  Interception  Scenario  Involving  Prohib¬ 
ited  Polygonal-Shaped  Areas. 

In  determining  if  I  is  capable  of  intercepting  T,  the  fol¬ 
lowing  simplifying  assumptions  are  made2: 

1 .  The  interceptor  I  and  the  threat  T  are  considered  to  be 
point  objects  flying  in  the  airspace  above  a  spherical 
earth  with  centre  O  and  radius  Re ; 

2.  /  and  T  fly  level  at  the  same  constant  altitude  h  and  at 
the  constant  speeds  Vi  =  \  \  vj(t)  \  |  and  vt  =  ||  Vx{t)  || 
respectively.  Moreover,  the  initial  position  Pt(0)  and 
velocity  v?{ 0)  of  T  are  known,  but  only  the  initial 
position  P/( 0)  and  speed  of  /  are  known  (the  initial 
direction  of  I  is  to  be  determined  as  part  of  the  solution 
to  the  intercept  problem); 

3.  The  maximum  weapons  range  of  I  is  Rw  and  the  time 
to  earliest  weapons  release  for  I  is  the  first  instant  that 
I  comes  within  great  circle  distance  Rw  of  T; 

4.  T  flies  along  the  great  circle  path  (trajectory)  deter¬ 
mined  by  its  initial  position  and  velocity.  Also,  the 
position  of  the  intercept  (if  it  is  possible)  and  the  posi¬ 
tion  of  I  at  all  times  lie  at  a  great  circle  distance  strictly 
less  than  i i(Re  +  h) /2  units  of  Pr(0); 

5.  I  is  to  intercept  T  in  the  minimum  time  possible.  It  is 
noted  that  since  /  is  assumed  to  fly  at  constant  speed, 
this  means  equivalently  that  I  is  assumed  to  intercept 
T  along  the  shortest  path  possible;  and 

6.  With  the  convention  that  a  point  in  the  airspace  lies 
above  an  area  on  the  earth’s  surface  if  and  only  if 
its  projection  onto  the  earth’s  surface  in  the  direction 
of  O  lies  in  the  area,  it  assumed  that  each  prohibited 
area  may  be  modelled  as  the  block  of  airspace  which 
lies  above  the  interior  of  a  spherical  polygon  on  the 
earth’s  surface.  Furthermore,  the  prohibited  areas  are 
assumed  to  be  pairwise  disjoint. 


2While  the  assumptions  have  been  made  either  to  assist  in  formulat¬ 
ing  the  problem  or  to  make  it  mathematically  tractable,  it  is  remarked  that 
none  of  the  assumptions  is  unreasonable  and  so  they  are  not  so  restrictive 
as  to  produce  misleading  results.  The  calculations,  which  are  made  at  each 
time  instant,  are  intended  to  support  decision-making  by  defence  operators 
who  would  be  monitoring  the  situation  as  it  unfolds  and  so  would  be  cog¬ 
nisant  of  any  substantial  changes  to  the  situation  that  might  invalidate  the 
assumptions  or  necessitate  a  new  evaluation. 


2.1  Mathematical  Preliminaries 

Throughout  the  remainder  of  the  paper,  the  following  math¬ 
ematical  notation,  terminology  and  results  are  employed. 

Unless  otherwise  specified,  standard  Euclidean  coordi¬ 
nates,  also  referred  to  as  earth-centred  rotating  (ECR)  co¬ 
ordinates  [3,  p.  176],  are  used  to  represent  points,  vectors, 
planes  and  all  other  geometrical  entities.  Vectors  are  writ¬ 
ten  either  as  lower  case  letters  in  boldface  such  as  v  or  in 
the  form  PQ  if  they  specify  the  vector  pointing  in  the  di¬ 
rection  of  point  Q  from  point  P.  In  particular,  the  position 
vector  of  the  point  P  is  denoted  by  OP  and  is  occasionally 
identified  with  P  for  the  sake  of  convenience.  For  all  non¬ 
zero  vectors  v,  the  unique  unit  vector  in  the  direction  of  v 
will  be  denoted  by  v.  All  other  standard  vector  notation, 
such  as  1 1  v  1 1  for  the  length  of  v,  is  also  adhered  to. 

A  continuous  function  /  :  [a,  b\  — >  R3  is  called  a  path 
in  M3  from  /(a)  to  f(b)  [4,  p.  59].  Unless  specified  oth¬ 
erwise,  all  paths  in  this  paper  are  assumed  to  be  smooth  al¬ 
most  everywhere  and  so  have  a  well-defined  length.  Given 
two  non-diametrically  opposed  points  P  and  Q  on  a  sphere 
of  radius  R  and  centre  ().  the  great  circle  defined  by  P  and 
Q  is  the  intersection  of  the  sphere  with  the  plane  determined 
by  P,  Q  and  O.  For  two  such  points  P  and  Q,  the  shorter 
of  the  two  segments  on  the  great  circle  from  P  to  Q  is  the 

great  circle  path  from  P  to  Q  and  is  denoted  by  PQ.  It 
can  be  shown  that  the  great  circle  path  from  P  to  Q  is  also 
the  shortest  path  on  the  sphere’s  surface  from  P  to  Q.  The 

great  circle  distance  between  P  and  Q  is  denoted  by  |  PQ\ 
and  is  related  to  the  Euclidean  distance  between  P  and  Q 
according  to  the  equation 

\\PQ\\  =  2Parcsin(||  PQ  ||  /2P).  (1) 

Furthermore,  setting  n  =  ||  OP  ||,  t  =  ||  v  ||,  where  v  = 

(OP  x  OQ)  x  OP,  and  d  =  | \PQ\ |,  a  general  point  on  the 
great  circle  path  between  P  and  Q  may  be  represented  in 
parametric  form  as  P(s)  such  that 

OP(s)  =  Rcos(s/R)n  +  Rsin(s/R)t ,  (2) 

where  s  £  [0,d]  denotes  the  arc  length  from  P  to  P(s) 
along  the  great  circle  path  (see  [5,  Eq.  58]).  It  is  also  noted 
that  n  is  normal  to  the  sphere  at  P  and  that  t  is  tangent  to 
the  sphere  at  P  in  the  direction  of  Q.  A  path  on  a  sphere  is 
said  to  be  polygonal  if  and  only  if  it  consists  of  a  finite  num¬ 
ber  of  great  circle  paths  such  that  no  two  connected  great 
circle  paths  lie  on  a  common  great  circle.  The  endpoints  of 
the  great  circle  paths  other  than  the  points  f(a)  and  f(b') 
are  referred  to  as  the  inner  vertices  of  the  path.  A  spherical 
polygon  is  a  simple  closed  polygonal  path  on  a  sphere  (that 
is,  a  simple  polygonal  path  for  which  f(b)  =  /(a)). 

A  graph  Q  consists  of  a  finite  set  V  and  a  finite  multiset 
£  consisting  of  two-element  multisets  from  V  [6,  p.  197], 
The  elements  of  V  are  called  vertices  and  those  of  £  are 
called  edges.  If  an  edge  e  contains  a  vertex  v,  then  e  is 
said  to  be  incident  with  v  and  v  is  said  to  be  an  endpoint 
of  e.  The  vertices  contained  in  a  given  edge  are  said  to  be 
adjacent.  Any  edge  which  has  the  same  vertex  for  both  its 
endpoints  is  called  a  loop  and  any  two  vertices  are  called 
multiply  connected  if  they  are  adjacent  with  respect  to  at 


least  two  distinct  edges.  Any  graph  without  loops  and  mul¬ 
tiply  connected  vertices  is  called  simple.  If  the  elements 
of  £  are  ordered  pairs,  then  Q  is  said  to  be  directed ,  other¬ 
wise  the  graph  is  said  to  be  undirected.  If  each  edge  has  a 
real  number  (a  weight )  assigned  to  it,  then  Q  is  said  to  be  a 
weighted  graph.  Finally,  a  sequence  of  vertices  and  edges 
Vl,  {vi,v2},v2,  {v2,f3},U3,  •  •  ■  ,Vn-l,{vn-l,Vn},Vn  for 
which  all  the  vertices  are  distinct  is  called  a  Q-path  from 
Vl  to  vn  3. 

3  Unobstructed  Threat  Interception 

As  a  first  step  towards  solving  the  complete  threat  inter¬ 
ception  problem,  means  of  solving  the  problem  in  the  ab¬ 
sence  of  prohibited  areas  are  investigated.  Throughout  this 
section,  it  is  assumed  that  the  threat  is  not  initially  within 
weapons  range  of  the  interceptor. 

Given  the  assumptions  on  the  dynamic  behaviour  and 
attributes  of  the  threat  T  and  the  interceptor  /,  the  unob¬ 
structed  threat  interception  problem  is  straightforward  to 
formulate.  Since  intercepting  T,  if  it  is  possible,  is  inter¬ 
preted  as  bringing  /  within  distance  R\y  of  T  in  the  mini¬ 
mum  time  possible,  it  follows  that  at  the  time  t*  of  the  inter¬ 
cept  the  great  circle  distance  between  T  and  P/(0)  will  be 
equal  to  the  distance  travelled  by  /  plus  its  weapons  radius, 
that  is 

\\PI(0)PT(t*)\\  =  vit*  +  Rw.  (3) 

Thus,  to  solve  the  intercept  problem,  it  is  sufficient  to  de¬ 
termine  the  smallest  positive  number  t*  satisfying  Eq.  3. 
Once  t*  is  found,  the  final  position  of  T  can  be  recovered 
via  Eq.  2  as 


OPr{t*)  =  Rcos{vTt*  /  R)ni  +  Rsin(i>Tt*  /  R)ti  (4) 


where  R  =  Re  +  h,  n\  =  OPt{ 0)  and  t\  =  'Wt(O).  This 
in  turn  allows  the  initial  velocity  of  I  to  be  determined  as 

•uT(0)  =  tt/it,  (5) 

where 

U  =  (oH( 0)  x  OPr(t*))  x  0^(0),  (6) 

and  the  final  position  of  I  to  be  determined  as 

OPi(t*)  =  R  cos(v  p*  /R)n2  +  Rsm(vit*  /  R)t2  (7) 

where  n2  =  OPj( 0),  t2  =  #/( 0)  =  u  and  R  =  Re  +  h  as 
before. 

Unfortunately,  it  is  impossible  to  solve  Eq.  3  in  closed 
form  in  general.  However,  fixed-point  algorithms  such  as 
the  Newton-Raphson  method  [7,  p.  79]  may  be  readily  ap¬ 
plied  to  find  t*  recursively  or  determine  that  no  such  such 
t*  exists,  as  the  case  may  be. 

A  second  approach,  which  is  only  valid  over  distances 
for  which  constant  velocity  motion  is  a  good  approximation 
to  motion  along  a  great  circle  path  at  constant  speed,  is  to 

3 It  is  more  customary  to  refer  to  such  a  sequence  simply  as  a  path  from 
vi  to  vn-  However,  the  term  (7-path  is  used  herein  to  avoid  any  possible 
confusion  between  a  path  in  the  graphical  sense  and  a  path  in  R3 . 


assume  that  T  and  /  both  travel  at  constant  velocity.  Under 
this  assumption,  Eq.  3  may  be  replaced  by  the  equation 

||  P/(0)/V(t*)  \\=vTt*  +  Rw,  (8) 

where  T  is  now  taken  to  be  within  weapons  range  of  I 
whenever  it  is  within  Euclidean  distance  R\y  of  /.  Unlike 
Eq.  3,  Eq.  8  is  straightforward  to  solve  in  closed  form.  On 
expanding  and  rearranging  the  terms  in  Eq.  8,  a  quadratic 
equation  a(f*)2  +  bt*  +  c  is  obtained  with 

a  =  Vj  —  Vt  (9) 

b  =  2(RwVl +  vT  ■  (01^(0) -OP^(0)) 

c  =  (R2w-  ||  OK(0)-O^(0)  ||2). 

A  detailed  analysis  of  the  quadratic  yields  the  following  so¬ 
lutions  for  the  smallest  positive  value  of  t*,  if  it  exists: 

{—c/b  if  a  =  0  and  b  >  0, 

(_6  +  AV2)/2a  if  a  <  0,  6  >  0 

and  A  >  0,  or  if  a  >  0, 
is  undefined  otherwise, 

where  A  =  62  —  4ac  denotes  the  discriminant  of  the 
quadratic.  Whenever  t*  exists,  the  final  position  of  T,  the 
velocity  of  I  and  the  final  position  of  I  are  respectively: 


OPT{t *) 

=  OPT(f))  +  vTt\ 

(11) 

vi 

=  u/P/(0)PT(f*), 

(12) 

oK(f*) 

=  oK(  0)  +  v/i*. 

(13) 

Otherwise,  the  intercept  is  impossible  and  so  these  terms 
are  undefined. 

4  Shortest  Path  Problem 

In  a  range  of  motion  planning  problems,  it  is  often  of  in¬ 
terest  to  determine  the  shortest  path  that  allows  an  object 
to  move  from  a  given  source  or  point  of  origin  to  a  given 
destination,  whilst  avoiding  particular  obstacles  that  exist 
in  the  object’s  environment  [8,  9,  10].  For  example,  in  ro¬ 
botics,  a  common  problem  is  that  of  determining  how  to 
plan  the  shortest  path  for  a  robot  to  follow  that  avoids  colli¬ 
sions  with  the  objects  in  the  room,  factory  or  terrain  etcetera 
in  which  it  is  operating.  Different  assumptions  about  the 
nature  of  the  object,  the  models  of  the  obstacles  and  the 
properties  of  the  desired  shortest  path  have  led  to  the  de¬ 
velopment  of  a  variety  of  approaches  and  treatments  of  the 
problem  in  the  plane  (see  [8,  9,  10,  11]  and  [12]  for  exam¬ 
ple).  The  approach  of  relevance  for  this  paper,  given  the 
assumptions  that  the  threat  and  interceptor  are  point  targets 
and  that  the  prohibited  areas  are  based  on  (spherical)  poly¬ 
gons,  is  that  which  employs  the  notion  of  a  visibility  graph 
and  Dijkstra’s  algorithm  for  finding  the  minimum  weight 
(/-path  (min  (/-path)  in  a  weighted  graph  [8,  9], 

4.1  Visibility  with  Respect  to  a  Planar  (Flat 
Earth)  Geometry 

In  the  context  of  the  shortest  path  problem  in  the  plane,  the 
notion  of  visibility  arises  naturally  through  the  need  to  de¬ 


termine  if  two  points  A  and  B  are  connected  by  a  line  seg¬ 
ment  that  does  not  intersect  any  of  the  obstacles.  In  general, 
points  A  and  B  are  defined  to  be  visible  to  each  other  with 
respect  to  a  given  set  S  if  no  point  of  S  lies  on  the  line  seg¬ 
ment  AB  \  { A ,  B}.  To  facilitate  the  solution  of  the  shortest 
path  problem  in  the  plane  involving  polygonal-shaped  ob¬ 
stacles,  in  which  A  and  B  are  the  point  of  origin  and  the 
destination,  a  graph  known  as  the  visibility  graph  may  be 
defined.  The  vertex  set  of  the  visibility  graph  comprises  the 
two  points  A  and  B  and  the  vertices  of  the  polygonal  ob¬ 
stacles.  Two  points  of  the  visibility  graph  are  then  adjacent 
by  definition  if  and  only  if  they  are  visible  to  each  other 
with  respect  to  the  obstacles  and  the  other  vertices  when  re¬ 
garded  as  points  in  the  plane.  The  visibility  graph  may  be 
used  to  construct  a  weighted  visibility  graph  by  attaching  a 
weight  to  each  edge  equal  in  value  to  the  Euclidean  distance 
between  the  two  endpoints  (regarded  as  points  in  the  plane) 
of  the  edge.  It  will  be  seen  in  Section  4.4  that  this  notion  of 
visibility  can  be  straightforwardly  generalised  for  solving 
the  shortest  path  problem  on  a  sphere.  For  example,  the  vis¬ 
ibility  graph  based  on  the  interceptor’s  initial  position  and 
the  set  of  vertices  of  the  polygons  in  Fig.  1  is  illustrated  in 
Fig.  2. 


Figure  2:  The  Visibility  Graph  Based  on  the  Interceptor’s 
Initial  Position  and  the  Polygon  Vertices  for  the  Target  In¬ 
terception  Example  in  Fig.  1. 

The  significance  of  the  visibility  and  weighted  visibility 
graphs  to  the  shortest  path  problem  in  the  plane  is  evident 
through  the  following  theorem. 

Theorem  1:  [9,  p.  309]  Any  shortest  path  between  A  and  B 
among  a  set  of  disjoint  polygonal  obstacles  is  a  polygonal 
path  whose  inner  vertices  are  vertices  of  the  obstacles. 

Theorem  1  implies  that  the  shortest  path  between  A  and 
B  which  avoids  the  obstacles  consists  of  a  sequence  of  line 
segments  whose  endpoints  (other  than  A  and  B )  are  (i)  ver¬ 
tices  of  the  obstacles  and  (ii)  visible  to  each  other.  Hence, 
the  shortest  path  in  the  plane  between  A  and  B  corresponds 
to  the  min  (/-path  in  the  weighted  visibility  graph. 

A  number  of  algorithms  of  differing  computational  com¬ 
plexity  exist  for  constructing  a  visibility  graph.  The  naive 
approach  to  the  problem  is  to  cycle  through  each  pair  of 
vertices  and  determine  their  visibility  with  respect  to  the 
edges  of  the  obstacles.  If  the  number  of  vertices  is  n,  then 
this  approach  leads  to  an  algorithm  that  runs  in  0(n3)  time 
since  the  number  of  pairs  of  vertices  is  0(n 2)  and  the  num¬ 
ber  of  edges  is  0(n)  (see  [8,  p.  297]  and  [9,  p.  310]).  Less 


computationally  expensive  algorithms  exist  [13],  but  for  the 
purposes  of  the  paper,  the  naive  approach  is  sufficient.  For 
a  more  in  depth  treatment  of  visibility  in  the  plane,  refer  to 
Chapter  19  of  Ref.  [10,  pp.  829-876], 


4.2  Dijkstra’s  Algorithm 

Having  constructed  the  weighted  visibility  graph,  it  still  re¬ 
mains  to  determine  the  min  [/-path  from  the  given  vertex  A 
to  the  other  given  vertex  B  to  find  the  shortest  path  from  A 
to  B.  One  well-known  and  often  applied  algorithm  for  de¬ 
termining  the  min  [7-path  between  two  vertices  is  Dijkstra’s 
algorithm,  which  was  discovered  by  the  Dutch  computer 
scientist  Edsger  Dijkstra  in  1959  [8],  In  fact,  Dijkstra’s  al¬ 
gorithm  determines  the  min  Q -path  from  the  first  given  ver¬ 
tex  A  to  every  other  vertex  in  the  graph.  Variants  of  Dijk¬ 
stra’s  algorithm  exist  for  directed  and  undirected  graphs  as 
well  as  for  simple  and  multiply-connected  graphs.  How¬ 
ever,  the  variant  of  relevance  here  is  that  which  applies  to 
simple,  undirected  graphs  since  weighted  visibility  graphs 
fall  into  this  class.  The  description  of  Dijkstra’s  algorithm 
presented  below  has  been  sourced  from  Ref.  [14]. 

For  convenience,  let  the  vertex  set  of  the  weighted  graph 
Q  be  V  =  {1, . . . ,  n}  with  1  as  the  given  vertex,  and  for 
each  i,  j  in  V  set  W (i,  j)  equal  to  0  if  i  =  j,  to  the  weight 
of  the  edge  {i,  j}  if  i  and  j  are  adjacent  (in  which  case  they 
are  necessarily  distinct)  and  to  oo  otherwise.  Then: 

Dijkstra’s  Algorithm 

SetL  =  {1} 

For  i  =  1, . . . ,  n 

Set  D{i)  =  W(l,i) 

If  W(l,t)  =  oo,  set  P(i)  =  0 
Otherwise  set  P(i)  =  1 
End  for 

While  V  \  L  ^  0 

Choose  k  £  V  \  L  with  D(k )  as  small  as  possible 

Put  k  in  L 

For  each  j  6  V\L 

If  D(j)>D(k)  +  W(k,j) 

Replace  D(j)  by  D (k)  +  W(k,  j) 
Replace  P(j)  by  k 

End  for 
End  while 
End 

At  the  completion  of  the  algorithm,  the  quantity  D(j)  is  the 
length  of  the  min  Q -path  from  1  to  j  (where  D(j)  =  oo  in¬ 
dicates  that  no  path  exists  between  1  and  j).  To  recover  the 
min  [7-path  itself,  the  function  P  can  be  used.  If  P(j)  =  0, 
then  no  path  exists  between  1  and  j.  Otherwise,  the  se¬ 
quence 


j,P(j),P(P(j)),P(P(P(j))),...  (14) 


For  example,  consider  the  visibility  graph  (7  depicted  in 
Fig.  2.  Representing  its  adjacency  matrix  as 


1  1  1  0  0  0  0  1  0 

1110  0  1110 
1  1  1  1  0  0  0  0  0 
0  0  1  1  10  0  0  0 
0  0  0  1111  0  1 
0  10  0  11111 
0  10  0  11111 
1  1  0  0  0  11  11 
_  0  0  0  0  1  1  1  1  1 


(15) 


where  the  first  row  and  column  correspond  to  the  intercep¬ 
tor’s  initial  position,  the  weighted  graph  Q\y  may  be  repre¬ 
sented  by 


0  321.7  417.6 

321.7  0  305.9 

417.6  305.9  0 

oo  oo  316.3 


306.9 

292.7 

278.3 


447.5 

418.5 


447.5  418.5 

0  58.4 

58.4  0 

454.2  406.8 

629.0  572.9 


543.8 

278.3 


454.2 

406.8 

0 

476.4 


417.3 
629.0 
572.9 

476.4 
0 


(16) 


Applying  Dijkstra’s  algorithm  to  the  weighted  visibility 
graph  Qw  above  yields  the  distance  matrix  D  and  the  path 
matrix  P  respectively  given  by: 


j 

D(j) 

P(j) 

1 

0 

1 

2 

321.7 

1 

3 

417.6 

1 

4 

733.9 

3 

5 

1032.9 

7 

6 

628.7 

2 

7 

614.4 

2 

8 

543.8 

1 

9 

1020.2 

8 

(17) 


Hence,  for  example,  the  shortest  distance  between  vertex  1 
(the  interceptor)  and  vertex  5  is  D( 5)  =  1032.9  and  the  17- 
path  along  which  it  is  attained  may  be  recovered  in  reverse 
order  as  5,  P( 5)  =  7,  P(7)  =  2,  P( 2)  =  1,  that  is  the  Q- 
path  is  1,  2,  7,  5. 


4.3  Shortest  Path  Problem  on  a  Sphere 

To  solve  the  shortest  path  problem  on  a  sphere  of  radius  R 
and  centre  O,  in  which  the  object  in  question  is  restricted 
to  move  along  a  piecewise  smooth  path  on  the  sphere,  the 
same  approach  used  for  the  planar  case  may  be  re-used 
with  minor  modifications  primarily  to  the  notion  of  visibil¬ 
ity.  Since  the  shortest  (piecewise)  path  between  two  (non- 
diametrically  opposed)  points  on  the  sphere  is  the  great  cir¬ 
cle  path  between  them  and  the  obstacles  are  spherical  poly¬ 
gons,  two  (non-diametrically  opposed)  points  A  and  B  on 
the  sphere  are  said  to  be  visible  to  each  other  with  respect 
to  a  set  S  if  and  only  if  no  point  of  S  lies  on  the  great  circle 

path  AB  \ { A ,  B}.  The  concepts  of  visibility  graphs  and 
weighted  visibility  graphs  are  then  defined  exactly  as  for  the 
planar  case,  and  so  Dijkstra’s  algorithm  also  remains  valid 
for  determining  the  min  Q -path.  However,  what  underpins 
the  generalisation  of  the  planar  approach  to  the  spherical 
case  is  the  fact  that  Theorem  1  remains  valid  if  “polygonal” 
is  replaced  with  “spherical  polygonal”,  that  is  the  following 
theorem  holds. 


lists  the  vertices  on  a  min  Q -path  from  1  to  j  in  reverse  Theorem  2:  Any  shortest  path  between  A  and  B  among  a 
order.  set  of  disjoint  spherical  polygonal  obstacles  is  a  spherical 


polygonal  path  whose  inner  vertices  are  vertices  of  the  ob¬ 
stacles. 

A  formal  proof  of  Theorem  2  established  by  the  authors 
will  appear  elsewhere.  However,  it  follows  the  same  logic 
as  that  employed  for  the  planar  case  (see  [8,  pp.  295-297] 
and  [9,  p.  309]).  Informally,  the  two  steps  in  the  proof  may 
be  summarised  as  follows.  First  it  can  be  shown  that  a  short¬ 
est  path  cannot  contain  a  smooth  subpath  other  than  a  great 
circle  path,  because  such  a  path  can  always  be  shortened 
by  replacing  part  of  it  with  a  great  circle  path.  Hence,  the 
shortest  path  consists  of  a  sequence  of  great  circle  paths. 
Second  it  can  be  shown  that  any  inner  vertex  of  the  path 
must  be  a  vertex  of  one  of  the  obstacles,  because  in  the 
neighbourhood  of  any  inner  vertex  that  were  not  also  a  ver¬ 
tex  of  one  of  the  obstacles,  the  path  could  be  shortened  by 
replacing  a  portion  of  it  by  a  suitable  great  circle  path.  As 
for  the  planar  case.  Theorem  2  guarantees  that  the  shortest 
path  corresponds  to  the  min  (/-path  in  the  weighted  visibil¬ 
ity  graph  because  consecutive  inner  vertices  are  visible  to 
each  other  (this  follows  from  the  facts  that  the  shortest  path 
by  definition  intersects  none  of  the  obstacles  and  the  seg¬ 
ment  of  the  shortest  path  between  them  has  been  shown  to 
be  a  great  circle  path). 

4.4  Visibility  with  Respect  to  a  Spherical 
Earth  Geometry 

Since  it  has  been  necessary  to  redefine  visibility  for  the 
spherical  case,  it  remains  to  demonstrate  means  of  deter¬ 
mining  the  visibility  graph.  To  exploit  the  existing  al¬ 
gorithms  for  the  planar  approach,  the  approach  proposed 
herein  is  to  project  each  point  of  the  set  S,  which  com¬ 
prises  the  point  of  origin,  the  destination  and  all  the  vertices 
of  the  obstacles,  onto  a  selected  tangent  plane  of  the  sphere 
in  such  a  way  that  visibility  is  preserved.  After  introduc¬ 
ing  local  coordinates  into  the  plane,  existing  algorithms  for 
the  planar  case  may  be  employed  to  construct  the  visibility 
graph.  The  details  follow. 

Projection 

To  systematically  choose  a  tangent  plane  into  which  to 
project  the  points,  all  that  is  needed  is  a  suitable  choice 
for  the  point  of  tangency.  One  solution  is  to  calculate  the 
barycentre  Pb  of  the  vertices  of  the  obstacles  and  use  its  pro¬ 
jection  onto  the  sphere  in  the  direction  of  the  vector  OPb 
as  the  point  of  tangency  r.  Given  m  obstacles  such  that  the 
ith  obstacle  has  rij  vertices  v,;;,  the  barycentre  Pb  is  given 
by 

m  rii  m 

Pb  =  EEwEn*  (lg) 

i—1  j—1  i—1 

and  so  the  point  of  tangency  r  is  given  by 

Or  =  ROPb.  (19) 

Then  each  point  P  of  S  is  projected  onto  the  plane  in  the 
direction  of  the  vector  OP.  Given  that  the  equation  of  the 
tangent  plane  at  r  is  given  by 

(20) 


a  simple  calculation  yields  the  projection  P'  of  P  onto  the 
plane  as 

OP' =  {R2/{Ch- -OP))OP.  (21) 

It  is  noted  that  the  projection  used  maps  the  great  circle 
paths  between  the  points  of  S  to  the  corresponding  line  seg¬ 
ments  between  the  projections  of  the  points  onto  the  tangent 
plane.  Therefore  two  points  in  S  are  visible  to  each  other  on 
the  sphere  if  and  only  if  their  projections  are  visible  to  each 
other  in  the  tangent  plane.  In  other  words,  the  projection 
preserves  visibility. 

Local  Coordinates 

Most  algorithms  for  determining  visibility  in  the  plane  rely 
on  algebraic  manipulations  of  the  2-dimensional  coordi¬ 
nates  of  the  vertices.  Therefore,  while  it  is  not  essential,  it 
is  arguably  preferable  to  transform  the  coordinates  of  points 
in  the  tangent  plane  from  their  3-dimensional  form  in  ECR 
coordinates  to  an  essentially  2-dimensional  form  in  local 
coordinates.  If  r  =  (0, 0,  ±P),  then  the  tangent  plane  has 
equation  z  =  ±P,  so  the  points  are  already  in  a  suitable 
form,  since  the  z  coordinate  of  each  point  is  the  same  and 
so  may  be  ignored.  Otherwise,  denote  r  by  (£1;  f2, £3)  and 
define  the  set  T  =  {ux,uy,uz}  where 


U Z  (*1,*2,*3) 

(22) 

Ux  =  (-*2,  *1,0) 

(23) 

Uy  =  Uz  X  Ux  • 

(24) 

Since  uz  (0,  0,  ±1),  it  is  straightforward  to  check  that 
the  vectors  in  T  form  an  orthonormal  basis  which  can  be 
used  to  define  the  axes  for  a  local  coordinate  system  with 
t  as  the  local  origin.  Since  uz  is  orthogonal  to  the  tan¬ 
gent  plane,  it  follows  that  ux  and  uy  lie  in  the  plane.  For 
any  point  x  =  (xi,X2,xf)  in  the  tangent  plane,  its  coor¬ 
dinates  with  respect  to  the  local  coordinates4  are  [x]r  = 
(xi,x'2,x'3)  where  (x  —  Or)  =  x[ux  +  x'2uy  +  x3uz. 
Taking  dot  products  of  both  sides  of  this  last  equation  with 
respect  to  each  of  the  basis  vectors  in  T  gives 

x\  =  (x  —  Or)  ■  ux  (25) 

x2  =  (x  —  Or)  '  uv  (26) 

x'3  =  (x  —  Or )  •  uz  =  0.  (27) 

Since  the  local  “z”  coordinate  is  0  for  each  point,  it  may 
be  ignored,  leaving  an  essentially  2-dimensional  coordinate 
system  as  required. 

Determining  Visibility 

Having  reduced  the  visibility  problem  on  the  sphere  to  the 
visibility  problem  in  the  plane,  any  existing  algorithm  for 
determining  visibility  may  be  used  to  construct  the  visibil¬ 
ity  graph.  For  discussions  on  how  to  construct  the  visibility 
graph  in  the  planar  case,  refer  to  [8]  and  [9]. 


4These  coordinates  are  sometimes  referred  to  as  East-North-Up  coor- 
dinates  because  of  the  directions  that  the  vectors  ux ,  uy  and  uz  point 
with  respect  to  a  spherical  earth  (cf.  the  North-East-Down  coordinates 
described  in  Ref.  [3,  p.  176]). 


(x  —  Or)  ■  Or  =  0 


5  Threat  Interception  in  the  Presence 
of  Prohibited  Areas 

The  techniques  and  results  from  Sections  3  and  4  may 
now  be  combined  to  determine  “satisficing”  solutions  to  the 
moving  threat  interception  problem  in  an  environment  con¬ 
taining  prohibited  areas. 

The  key  to  the  solution  is  the  observation  that  if  the  in¬ 
terceptor  /  is  able  to  intercept  the  threat  T  and  that  it  does 
so  in  the  shortest  time  possible,  then  the  final  great  circle 
path  (leading  to  the  weapon  release  point)  emanates  either 
from  the  initial  position  of  I  or  from  a  vertex  of  one  of 
the  prohibited  areas.  Thus,  the  shortest  path  to  intercept 
T,  if  one  exists,  lies  in  the  set  Sn  of  paths  constructed  as 
follows.  Let  S  denote  the  set  consisting  of  the  initial  posi¬ 
tion  P/(0)  of  /  and  the  n  vertices  of  the  prohibited  areas. 
For  each  element  V  of  S,  compute  the  time  ty  required 
by  I  to  travel  along  the  shortest  path  7Tmin(U)  from  Pj( 0) 
to  V  and  determine  if  /  can  intercept  T  (now  at  position 
Pr(tv))  from  V  along  a  path  such  that  T  remains  visi¬ 
ble  to  I.  If  so,  for  each  such  path,  define  the  augmented 

path  7Tmin (Vj  U  VPj  ( ty ),  where  ty  is  the  time  to  the 
intercept  from  V.  As  V  varies  over  the  set  S,  the  union 

of  all  such  paths  Try  =  7Tmin(U)  U  VPi  (ty)  constitutes 
the  set  ,SL .  Then  I  is  able  to  intercept  T  whilst  travelling 
at  constant  speed  ||  vj  ||  and  also  avoiding  the  prohibited 
areas  if  and  only  if  .S'-  is  non-empty.  The  heuristic  em¬ 
ployed  here  to  find  reasonable  (suboptimal)  solutions  is  not 
to  consider  every  path  in  Sn,  but  to  restrict  attention  to  those 
paths  which  are  constructed  as  outlined  in  the  following  al¬ 
gorithm.  From  amongst  those  paths  considered,  the  path 
ultimately  selected  (assuming  at  least  one  exists)  is  the  one 
for  which  the  time  to  earliest  weapon  release  is  minimum. 

Initialisation  Set  the  minimum  time  to  intercept  (that  is, 
the  time  to  earliest  weapon  release)  to  oo  and  the  shortest 
intercept  path  from  Pj( 0)  to  T  to  0. 

Step  1:  Given  the  initial  conditions  for  the  threat  T  and  the 
interceptor  /,  determine  if  I  can  intercept  T  assuming  that 
the  prohibited  areas  do  not  exist.  If  it  cannot,  then  halt  the 
algorithm5.  Else  proceed  to  step  2. 

Step  2:  Still  under  the  same  assumption,  check  if  the  result¬ 
ing  path  that  I  would  follow  to  intercept  T  passes  through 
any  of  the  prohibited  areas  (that  is,  determine  the  visibility 
of  F s  initial  position  and  its  final  position  with  respect  to 
the  prohibited  areas).  If  the  initial  and  final  positions  are 
visible  to  each  other  with  respect  to  the  polygons,  then  set 
the  minimum  time  to  intercept  to  the  time  required  for  I  to 
travel  along  this  great  circle  path,  set  the  shortest  intercept 
path  to  the  great  circle  path6  and  halt  the  algorithm.  Other¬ 
wise,  proceed  to  step  3. 

Step  3:  Using  the  techniques  described  in  Section  4,  con- 

5lf/  cannot  intercept  T  under  this  assumption,  then  it  will  not  be  able 
to  intercept  T  in  the  presence  of  the  prohibited  areas  either  because  any 
path  that  I  were  to  take  around  the  prohibited  areas  would  necessarily  be 
longer  than  the  great  circle  path  between  T  and  I. 

6  In  this  case,  I  can  intercept  T  along  a  direct  great  circle  path,  so  that 
the  prohibited  areas  provide  no  obstacle.  Since  any  other  path  around  the 
prohibited  areas  would  necessarily  be  longer  than  the  great  circle  path,  the 
shortest  intercept  path  in  this  case  has  already  been  found. 


struct  the  visibility  graph  Q  for  the  set  of  vertices  S  consist¬ 
ing  of  Pj( 0)  and  the  n  vertices  of  the  prohibited  areas7. 

Step  4:  Construct  the  weighted  visibility  graph  Gw  by  as¬ 
sociating  a  weight  to  each  edge  of  Q  that  is  equal  in  value  to 
the  great  circle  distance  between  the  corresponding  vertices 
of  S.  Then  apply  Dijkstra’s  algorithm  to  Gw  to  calculate 
the  min  (/yp-paths  from  P/(0)  to  each  of  the  other  vertices 
and  the  lengths  of  these  paths. 

Step  5:  For  each  vertex  V  €  S  \  {P/(0)},  if  the  shortest 
path  from  Pj( 0)  to  V  is  finite: 

1.  Determine  the  time  ty  required  for  I  to  travel  to  V 
along  the  path  and  determine  the  position  of  T  at  ty. 

2.  Determine  if  /  at  position  V  can  intercept  T  (now  at 
position  Pr(ty))  at  some  time  ty,  under  the  assump¬ 
tion  that  there  are  no  prohibited  areas.  If  so,  determine 
if  Pr(ty)  lies  within  a  great  circle  distance  nR/2  of 
Pt( 0),  and  that  V  and  Pr(ty)  are  visible  to  each 
other.  If  so,  then  if  (ty  +  ty)  is  less  than  the  cur¬ 
rent  minimum  time  to  intercept,  reset  the  minimum 
time  to  intercept  to  (ty  +  ty)  and  the  shortest  path 

to  TTy  =  7Tmin(U)  U  UP/  (ty). 

Once  all  the  steps  have  been  completed,  if  the  minimum 
time  to  intercept  is  infinite,  then  I  is  not  able  to  intercept  T 
using  the  adopted  strategy.  Otherwise,  I  is  able  to  intercept 
T  and  the  algorithm  returns  an  efficient  path  from  P/  (0)  to 
the  point  of  intercept  and  the  time  to  earliest  weapon  release 
along  this  path. 

6  Example 

Consider  the  interception  scenario  illustrated  in  Fig.  1 .  The 
interceptor’s  speed  was  1000  km/hr  and  it  had  a  weapons 
radius  of  30  km.  The  threat’s  speed  was  700  km/hr.  At  step 
1,  using  the  first  unobstructed  threat  interception  technique 
described  in  Section  3  resulted  in  a  potential  interception 
solution.  However,  at  step  2,  the  initial  and  final  interceptor 
positions  were  found  not  to  be  visible  to  each  other  with 
respect  to  the  polygonal  prohibited  areas,  and  so  the  algo¬ 
rithm  proceeded  to  steps  3  and  4.  At  these  steps,  the  visi¬ 
bility  graph  Q  and  the  weighted  visibility  graph  Gw  in  Eqs. 
15  and  16  were  constructed,  and  Dijkstra’s  algorithm  was 
employed  to  find  the  shortest  distances  and  the  associated 
Gw- paths  in  Eq.  17  from  the  initial  interceptor  position 
to  each  vertex  of  the  prohibited  (polygonal)  areas.  Finally, 
at  step  5,  exactly  two  threat  interception  solutions  resulted. 
The  first  path  proposed  was  1,  2,  7,  5  to  the  threat  and  the 
second  path  proposed  was  1,  8,  9  to  the  threat.  The  first 
path  was  the  shorter  of  the  two,  so  the  output  of  the  algo¬ 
rithm  was  the  path  1,  2,  7,  5  to  the  threat  indicated  in  green 
in  Fig.  3  with  a  time  to  earliest  weapon  release  of  8.43 
minutes  (the  figure  also  indicates  the  strike  envelope  of  the 

7Note  that  it  is  not  necessary  to  construct  the  whole  visibility  graph 
from  scratch  at  each  iteration.  Instead,  the  (visibility)  subgraph  of  Q  based 
on  just  the  vertices  of  the  prohibited  areas  can  be  constructed  once,  during 
the  initialisation  stage  for  example.  Then,  at  each  iteration,  it  is  only  nec¬ 
essary  to  determine  the  visibility  of  Pj  (0)  to  each  of  the  vertices  of  the 
prohibited  areas  to  complete  the  subgraph  to  Q.  This  substantially  reduces 
the  computational  expense  of  constructing  Q. 


interceptor  by  the  green  circle  and  the  threat’s  position  at 
the  time  of  the  intercept  by  the  red  asterisk). 


Figure  3:  The  Most  Efficient  Path  from  the  Interceptor  to 
the  Threat  as  Determined  by  the  Algorithm. 


7  Conclusion 

In  the  context  of  weapons  allocation  for  tactical  air  defence, 
the  problem  of  threat  interception  in  an  environment  con¬ 
taining  prohibited  areas  has  been  investigated  and  an  ap¬ 
proach  based  on  techniques  from  the  discipline  of  compu¬ 
tational  geometry  has  been  proposed.  In  particular,  by  inter¬ 
preting  threat  interception  as  a  shortest  path  problem,  solu¬ 
tions  have  been  developed  for  intercepting  static  and  mov¬ 
ing  threats  by  extending  a  known  approach  to  the  shortest 
path  problem  in  the  plane,  based  on  visibility  graphs  and 
Dijkstra’s  algorithm.  For  static  threats,  the  solution  based 
on  the  generalisation  of  the  notion  of  visibility  for  a  spher¬ 
ical  earth  geometry  yields  the  shortest  path  from  the  inter¬ 
ceptor  to  the  threat,  while  for  moving  threats,  the  approach 
described  is  “satisficing”  in  nature,  producing  suboptimal, 
yet  efficient  solutions.  An  example  has  also  been  given  to 
illustrate  the  ideas  developed  herein.  Finally,  it  is  noted 
that  while  this  problem  has  been  discussed  in  isolation  of 
the  threat  evaluation  process  which  led  to  determining  that 
the  threat  should  be  intercepted,  the  two  processes  are  bet¬ 
ter  regarded  as  different  interdependent  aspects  of  a  single 
threat  evaluation  and  weapons  allocation  (TEWA)  process 
[15,  16].  As  such,  future  work  will  focus  on  incorporating 
this  threat  interception  technique  into  a  TEWA  testbed  for 
the  purpose  of  experimenting  with  emerging  TEWA  con¬ 
cepts,  as  well  as  testing  and  evaluating  different  threat  in¬ 
terception  strategies. 
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